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Previous  studies  have  indicated  that  a  novel  design  for  proton  exchange  membrane  fuel  cells  with  a 
piezoelectric  (PZT)  device,  which  is  regarded  as  an  actuator  for  pumping  air  onto  the  cathode  channel,  can 
offer  better  performance  with  higher  current  generation.  These  results  indicate  that  piezoelectric  proton 
exchange  membrane  fuel  cells  (PZT-PEMFCs)  may  compress  more  air  into  the  catalyst  layer  and  thus  may 
enhance  electrochemical  reactions,  resulting  in  higher  current  output.  At  the  same  time,  produced  water 
vapor  is  pumped  out  from  the  cathode  channel  during  the  compression  process.  Previous  studies  on  PZT- 
PEMFCs  without  ribs  showed  the  strong  effect  of  ohmic  and  concentration  losses.  In  this  study,  a  shallow 
rib  is  chosen  to  reduce  the  aforementioned  losses  and  pressure  drop  in  the  cathode  channel.  The  rib  design 
is  an  important  parameter  that  can  be  used  as  the  support  for  the  membrane  electrode  assembly  (MEA).  A 
transient  three-dimensional  model  is  built  to  simulate  and  compare  the  performance  of  PZT-PEMFCs  both 
with  and  without  ribs.  Water  vapor,  oxygen,  and  current  density  profiles  in  the  PZT-PEMFC  are  studied  in 
detail.  The  major  operating  parameters  include  the  rib  width  and  the  PZT  vibration  frequency.  Our  results 
show  that  the  ribbed  cathode  channel  can  reduce  ohmic  losses  and  double  current  generation.  Moreover, 
at  higher  PZT  vibration  frequency  (f=  64  Hz),  an  air-breathing  PZT-PEMFC  compresses  more  oxygen  into 
the  catalyst  layer  and  thus  enhances  the  electrochemical  reaction,  resulting  in  a  higher  current  output 
(0.208  A  cm-2). 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Air-breathing  proton  exchange  membrane  fuel  cells  (AB- 
PEMFCs)  have  attracted  attention  as  substitutes  for  currently  used 
batteries.  During  AB-PEMFC  operation,  liquid  water  is  transported 
by  electro-osmotic  drag,  back  diffusion,  and  convection  in  the 
membrane.  The  electro-osmotic  drag  coefficient  has  been  found  to 
be  a  function  of  the  water  content  of  the  membrane  and  tempera¬ 
ture  [1-4].  Generally,  the  electro-osmotic  drag  coefficient  increases 
with  increasing  water  content  in  the  membrane.  At  high  water  con¬ 
tent,  the  electro-osmotic  drag  coefficient  increases  with  increasing 
temperature.  Yi  and  Nguyen  [5]  showed  that  PEMFC  performance 
was  improved  by  anode  humidification  and  positive  differential 
pressure  between  the  cathode  and  anode.  Ge  and  Yi  [6]  predicted 
that  the  dry  reactant  gases  could  be  successfully  humidified  inter¬ 
nally  and  could  maintain  high  performance  when  the  PEMFC  was 
operated  in  counter  flow  mode.  Thus,  water  transport  and  man¬ 
agement  is  dependent  upon  the  structure  and  properties  of  the  cell 
components,  reactant  stream  humidification,  flow  field  layout,  and 
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structural  and  wetting  properties  of  the  gas  diffusion  media  and 
micro-porous  layer  [7]. 

In  general,  AB-PEMFCs  do  not  need  an  additional  oxygen  supply 
device;  their  performance  is  proportional  to  the  oxygen  transfer 
coefficient  [8].  Santa  Rosa  et  al.  [9]  developed  an  open-air  cathode 
PEMFC  stack  with  an  air  fan.  This  design  combined  oxidant  sup¬ 
ply  and  stack  cooling.  Using  numerical  analysis,  Zhang  et  al.  [10] 
presented  various  stack  and  cell  parameters  for  stack  designs  with 
an  array  of  air-breathing  fuel  cell  cartridges.  Matamoros  and  Brug- 
gmann  [11  ]  showed  that  dehydrating  phenomena  slightly  affected 
the  performance  of  AB-PEMFCs.  Moreover,  catalyst  loading,  rela¬ 
tive  humidity,  temperature,  hydrogen  stoichiometry,  gas-diffusion 
layer  thickness,  and  cathode  structure  are  important  parameters  in 
the  performance  of  AB-PEMFCs  [12,13]. 

Proper  flow  field  design  of  the  AB-PEMFC  may  help  address  the 
water-flooding  problem  and  enhance  cell  performance.  For  inter- 
digitated  flow  fields,  the  reactant  gas  can  be  delivered  to  the  catalyst 
layer  and  can  increase  the  electrochemical  reaction  rate  [14].  How¬ 
ever,  while  pin-type  flow  fields  have  a  lower  reactant  pressure, 
they  have  uneven  gas  distribution  and  poor  water  removal  and  cell 
performance  [15].  The  single-serpentine  flow  field  exhibits  better 
performance  than  double-  and  triple-serpentine  flow  fields  [16]. 
In  addition,  a  new  convection-enhanced  serpentine  flow  field  has 
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Nomenclature 

Ain  inlet  area  (m2) 

APZT  piezoelectric  area  (m2 ) 

Deff  effective  diffusion  coefficient  (m  s-1 ) 

/  frequency  of  PZT  (Hz) 

F  Faraday  constant  (96,439  C  mole-1 ) 

h  enthalpy  (kj) 

i  current  (A) 

Ji  species  diffusive  flux  (mol  m-2  s-1 ) 

jx  transfer  current  (A) 

jo  reference  current  density  at  a  known  open  circuit 

voltage  (A  m-3 ) 

ke ff  effective  thermal  conductivity  (W  m-1  K-1 ) 

M  mass  (kg) 

N  amount  of  species 

P  pressure  (Nm-2) 

R  gas  constant  (J  mol-1  K-1 ) 

t  time  (s) 

T  temperature  (K) 

Va  anode  inlet  velocity  (m  s-1 ) 

Vc  cathode  inlet  velocity  (m  s-1 ) 

\/pzT  motion  equation  of  the  piezoelectric  device  (m  s-1 ) 
Yt  mass  fraction  of  ith  species 

v  volume  (m3) 

Greek  letters 
s  porosity 

r\  over  potential  (V) 

k  permeability  (m2) 

fi  viscosity  (kg  m-1  s-1 ) 

p  density  (kg  m-3) 

r  shear  stress  (N  m-2 ) 

0  equivalence  ratio 

oil  production  rates  (kg  m-3  s-1 ) 

Subscript 

act.  actual 

stoi.  stoichiometry 


been  proposed  by  Xu  and  Zhao  [17]  to  increase  the  mass  transport 
rate  and  remove  liquid  water  trapped  in  the  porous  media.  Further¬ 
more,  the  parallel  flow  field  with  an  incorporated  baffle  exhibited 
better  performance  because  the  baffle  effect  forces  the  reactant 
gas  to  pass  through  the  diffusion  layer  [18].  Moreover,  Ma  et  al. 
[19,20]  have  proposed  a  novel  cathode  channel  design  that  utilizes 
the  piezoelectric  (PZT)  effect  in  AB-PEMFC  systems,  known  as  the 
PZT-PEMFC.  The  PZT-PEMFC  can  feed  air  into  the  cathode  chan¬ 
nel  through  a  gas  pump  with  a  piezoelectric  actuation  structure.  A 
higher  PZT  vibration  frequency  leads  to  a  higher  oxygen  concen¬ 
tration  at  the  cathode  channel,  and  most  of  the  produced  water 
vapor  can  be  pumped  out  of  the  cathode  channel  by  PZT  vibrations. 
This  design  can  reduce  the  water-flooding  problem  and  increase 
cell  efficiency. 

Studies  have  shown  that  a  ribbed  design  affects  cell  perfor¬ 
mance  and  species  distribution,  and  that  a  narrower  channel  with 
a  wider  rib  spacing  yields  better  performance  [21,22].  Although 
PZT-PEMFCs  may  increase  cell  efficiency  and  reduce  water  flood¬ 
ing,  a  rib-free  design  may  induce  a  large  resistance.  Thus,  ribs 
are  necessary  on  the  cathode  side  of  the  PZT-PEMFC.  Jeong  et  al. 
[23]  showed  that  a  single  cell  with  a  cathode  open  area  of  77% 
exhibited  the  best  performance.  This  result  was  related  to  oxy¬ 
gen  transport  from  the  atmosphere  to  the  catalyst  layer.  Hottinen 


et  al.  [24]  concluded  that  the  resistance  of  a  larger  open  cath¬ 
ode  area  increased  tremendously,  adversely  affecting  performance. 
Thus,  a  novel  cathode  side  design  with  ribs  in  the  PZT-PEMFC 
is  proposed  in  this  study.  In  addition,  a  valveless  design  with 
four  air  ports  is  applied  in  the  PZT-PEMFC,  differing  from  previ¬ 
ous  studies  [19,20].  A  transitional  three-dimensional  model  [25], 
based  on  the  semi-implicit  method  for  pressure  line  equations 
consistent  (SIMPLEC)  procedure,  is  used  to  calculate  the  effect  of 
ribs  on  the  performance  of  the  PEMFC  at  different  PZT  frequen¬ 
cies. 

2.  Mechanisms  of  PZT-PEMFC 

The  PZT-PEMFC  (Fig.  1)  contains  a  piezoelectrically  actuated 
micro-diaphragm  that  is  installed  in  the  cathode  channel.  This 
device  is  regarded  as  an  original  gas  pump  design  and  can  draw 
air  into  the  cathode  channel  without  an  external  air  pump.  While 
the  actuator  moves  outward,  the  volume  of  the  cathode  channel 
increases,  and  the  air  is  drawn  into  the  chamber  (Fig.  2).  How¬ 
ever,  when  the  actuator  moves  inward,  the  volume  of  the  cathode 
channel  decreases,  and  the  air  is  compressed  into  the  catalyst  layer. 
Previous  studies  have  proposed  an  ideal  air-breathing  PZT-PEMFC 
cycle  that  describes  four  piezoelectric  actuation  processes.  In  addi¬ 
tion,  the  motion  of  the  PZT  device  was  assumed  to  be  a  step  function 
[19,20].  In  this  study,  a  sine  function  was  applied  for  the  motion  of 
PZT  device  as  it  may  improve  both  the  electrochemical  reaction  in 
the  catalyst  layer  and  the  power  output. 

In  addition,  the  PZT-PEMFC  can  also  be  used  in  stacks.  Fig.  3 
shows  that  two  cathode  channels  and  two-cell  stacks  can  be  driven 
by  one  PZT  device. 

2.1.  Equivalence  ratio 

The  equivalence  ratio,  defined  as  the  actual  fuel-to-air  ratio 
divided  by  theoretical  fuel-to-air  ratio,  is  an  important  parameter 


Air 

Hydrogen  Inlet 


PZT  Device 


Air 

Hydrogen  Outlet 


Fig.  1.  The  outlook  of  a  PZT-PEMFC  system. 
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MEA 


(a)  Sucked  air  into  the  channel 
MEA 


Fig.  2.  Actuating  mechanisms  in  a  PZT-PEMFC  system. 


Fig.  4.  The  sine  function  of  periodical  flow  at  cathode. 


the  cathode  is  expressed  by  the  Reynolds  Transport  Theorem  as 

DM  _  D(pV) 

~Dt  ~  Dt 

=  |-  /  P  dV  +  /  pV m  It  dA+  pc\tc  ■  ~ft  dA  =  0. 

^  Jcv  JCS,  PZT  Jcs,  in 

(2) 


for  discussing  air-fuel  mixing  phenomena  in  internal  combustion 
engines  [26].  Similarly,  the  equivalence  ratio  is  also  an  important 
parameter  for  air-breathing  PZT-PEMFCs.  The  theoretical  demand 
of  air  for  generating  1A  of  current  is  17.96  mg  min-1,  and  the 
(hydrogen/air)stoi  is  0.03.  In  this  study,  the  equivalence  ratio  0 
is  defined  as  the  ratio  of  ( hydrogen/air )act.  to  ( hydrogen/air )stoi 
Fuel-rich  conditions  (0  >  1 )  are  found  in  air-deficient  and  incom¬ 
plete  hydrogen  reactions  at  low  PZT  vibration  frequencies.  On 
the  contrary,  fuel  lean  conditions  (0<1)  are  found  at  high  PZT 
vibration  frequencies,  where  sufficient  air  is  drawn  into  the 
cathode  channel,  resulting  in  an  increase  in  hydrogen  consump¬ 
tion. 


The  gas  density  is  expressed  by  the  equation  of  state  for  an  ideal 
gas,  so  the  mass  conservation  equation  is  derived  as 


Rdt  T  ^  +  Pcs,pztWzt^pzt  -  PcVcA n  =  0. 

J  cv 

The  inlet  velocity  at  the  cathode,  V^,  is  derived  as 


v£  = 


1 

PcAn 


11  [  P 
RdtJc,T 


dV  +  pcs VpzjApzr 


(3) 


(4) 


3.  Theoretical  model  for  PZT-PEMFC 


2.2.  Actuation  analysis  of  PZT-PEMFC 

The  driving  forces  at  the  anode  and  cathode  inlets  are  different. 
The  inlet  hydrogen  velocity  at  the  anode  is  assumed  to  be  constant. 
However,  the  inlet  air  velocity  at  the  cathode  is  driven  by  PZT  vibra¬ 
tions.  Thus,  the  velocity  V  is  divided  into  the  anode  velocity  Va  and 
cathode  velocity  Vc.  The  equation  of  motion  for  the  PZT  is  a  sine 
function  given  by 

^  =  s{“a0005  x  (sin(2?r/t  ■?))}•  (1) 

Also,  the  inflow  and  outflow  periods  in  the  channel  as  induced 
by  the  sine  function  are  shown  in  Fig.  4.  The  inlet  air-flow  rate  at 


Anode  Channel 
MEA 

Cathode  Channel 

PZT  Dev  ice 
Cathode  Channel 
MEA 

Anode  Channel 
Fig.  3.  The  outlook  of  a  PZT-PEMFC  stack. 


A  numerical  model  of  the  novel  fuel  cell,  based  on  the  SIM- 
PLEC  procedure,  has  also  been  developed.  The  assumptions  are  as 
follows: 

(1)  The  reactants  and  products  are  treated  as  ideal  gases. 

(2)  The  Stefan-Maxwell  equations  are  applied  to  multispecies  dif¬ 
fusion. 

(3)  Contact  electrical  and  heat  resistances  among  the  channel 
layer,  diffusion  layer,  catalyst  layer,  and  membrane  are  ignored. 

(4)  The  Nernst-Planck  equation  is  used  for  proton  transport 
through  the  membrane. 

(5)  Ohm’s  law  is  applied  across  the  entire  region  of  the  polariza¬ 
tion  curve. 

(6)  The  porosity  and  permeability  of  the  porous  media  are  uni¬ 
form,  and  the  porous  media  are  isotropic  and  homogeneous. 

(7)  The  effect  of  gravity  is  ignored. 

(8)  The  amplitude  of  the  PZT  device  is  assumed  constant  for  all 
frequencies,  and  the  motion  of  PZT  device  is  assumed  to  be  a 
sine  function. 

(9)  The  water  in  the  fuel  cell  is  assumed  to  be  vapor. 

(10)  The  relative  humidity  of  fuel  and  air  is  100%. 

The  membrane  electrode  assembly  (MEA)  parameters  refer  to  a 
DuPont  membrane,  and  the  dimensions  of  the  PEMFC  are  shown 
in  Fig.  5.  The  compositions  of  the  inlet  humidified  fuel  and  air  are 
shown  in  Table  1 .  In  addition,  the  inlet  fuel  velocity  is  1  m  s-1 .  The 
inlet  air  velocities  vary  with  different  PZT  frequencies,  as  shown  in 
Table  2. 
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Fig.  5.  Dimensions  of  PZT-PEMFC. 


Table  1 

Gas  composition  of  fuel  and  air  at  inlet  conditions 


Fuel 

Hydrogen 

43.9% 

Water  vapor 

56.1% 

Nitrogen 

69.8% 

Air 

Oxygen 

22.1% 

Water  vapor 

8.1% 

Unit:  molar  concentration. 


3.2.  Diffusion  layer 

For  porous  media,  the  governing  equations  should  be  modified 
by  the  porosity  e  and  the  permeability  k  of  the  media.  The  governing 
equations  are  as  follows: 

The  Continuity  equation  is 

|(ep)  +  V.(epV)  =  0.  (10) 

The  momentum  equation  is 

fj  ^  p2uV 

^(spV)  +  (V  •  spVV)  =  -sVp  +  V  •  {ex)  +  — (11 ) 

where  £2/xV//<  is  Darcy’s  drag  force  describing  the  fluid  flowing 
through  the  porous  media. 

The  energy  equation  is  written  as 

^(eph)  +  V  •  (. epVah )  =  V  ■  q  +  ez  :  VVa  +  -jTr;  +  '-f.  (12) 

The  species  equation  is 

j|(«pYf)  +  V  •  {spVaYt)  =  V  Jj.  (13) 


3.1  Channel  layer 


3.3.  Catalyst  layer 


The  governing  equations  in  the  channel  include  the  continuity, 
momentum,  and  species  equations.  The  continuity  equation  is 

|?  +  V.(pV)  =  0.  (5) 

The  momentum  equation  is 

2(pV)  +  (W-pVV)  =  -Vp  +  W-r.  (6) 

The  energy  equation  is 

2(ph)  +  V-(pVh)  =  V-q  +  r:  VV+^.  (7) 


The  important  governing  equations  in  the  catalyst  layer  are  the 
species  and  current  conservation  equations.  Because  the  electro¬ 
chemical  reaction  and  current  generation  take  place  in  this  layer, 
the  species  and  current  conservation  equations  should  be  written 
in  the  forms  given  below. 

The  species  equation  is 

%(epY, )  +  V  ■  ( epVYt )  =  V  J,  +  oj„  (14) 

where  is  the  production  rate  in  the  gas  phase  of  the  ith  species. 
The  transfer  current  jT  can  be  expressed  by  the  Butler-Volmer  equa¬ 
tion  as 


The  species  equation  is 

+  V  •  (fiVYi)  =  V  (8) 

where  Yz  is  the  mass  fraction  in  the  gas  phase  of  the  ith  species.  Ji 
is  the  diffusive  flux,  which  can  be  expressed  as 

h  =  pD,effVy,  +  ^DwffVM  -  effVYj 

j 

-pVMj^DteffVYj,  (9) 

j 

where  pDieffVYi  is  the  Fickian  diffusion  coefficient,  which  is 
induced  by  the  concentration  gradient.  Other  terms  are  the  cor¬ 
rections  of  multiple  species  within  the  porous  medium.  The  fluid 
diffusion  coefficient  should  be  modified  to  Dz  eff,  which  is  derived 
by  Bruggeman  [27]. 


Table  2 

Operating  parameters 


1  ms-1 

Anode  inlet 

Pressure 

1  atm 

Temperature 

323  K 

Vc 

Eq.  (4),  m  s_1 

Cathode  inlet 

Pressure 

1  atm 

Temperature 

323  K 

J  T  = 


Jo _ 

nN  r  a  1 

ref] 


exp 


^nVexp^ 


ak 


V  RT 


V  RT 


0]D*i  • 

k=  1 


(15) 


where  p  is  an  over  potential,  defined  as  the  difference  between  the 
solid  potential  and  the  porous  (or  liquid)  potential. 


Fig.  6.  Current  density  profiles  with/without  rib  (f=  8  Hz). 
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t  =  0.5  sec. 


H20 

0561 


Fig.  7.  Water  vapor  distribution  in  different  time  (rib  width  =  2  mm, /=  8  Hz). 


Fig.  8.  Current  density  profiles  with  different  rib  width  (J=  16  Hz). 

The  current  conservation  equation  is  given  by 

V  i  =  0.  (16) 

Because  the  current  flows  through  the  porous  medium,  it  can  be 
divided  into  two  parts.  The  first  part  is 

I  =  if  +  Is,  (17) 

where  if  is  the  current  flowing  through  the  fluid,  and  is  is  the  cur¬ 
rent  flowing  through  the  solid.  Because  they  are  porous  media,  the 


electrons  transfer  from  fluid  to  solid  and  solid  to  fluid.  Thus,  the 
transfer  current  can  be  written  as 

> 

i 

II 

> 

II 

E— 1 

(18) 

Using  Ohm’s  law,  this  can  be  changed  to 

V  •  (crfV (fif)  =  -V  •  (crsV<ps)  =jT. 

(19) 

4.  Results  and  discussion 

A  novel  ribbed  PZT-PEMFCs  design  has  been  developed,  and  a 
three-dimensional,  transitional  model  has  been  successfully  built 
to  analyze  its  performance  and  primary  phenomena. 

4.1  Influence  of  rib  design  on  current  density 

The  rib  design  is  an  important  parameter  in  the  fuel  cell  for 
reducing  internal  resistance  and  for  preventing  MEA  fracture  due  to 
the  pressure  difference  between  the  anode  and  cathode.  However, 
previous  studies  [19,20]  did  not  have  any  ribs  in  the  cathode  chan¬ 
nel.  As  shown  in  Fig.  6,  only  the  unribbed  design  has  a  small  current 
density,  0.12  A  cm-2,  due  to  higher  resistance.  A  ribbed  cathode 
channel  can  reduce  cathode  open  air  area,  which  decreases  the 
resistance  and  doubles  the  current  density  [22,23].  However,  the 
ribbed  design  has  a  water-accumulation  problem.  Therefore,  the 
current  decreases  with  time.  As  shown  in  Fig.  7,  water  vapor  may 
be  blocked  and  may  be  accumulated  near  the  rib  with  time. 


Fig.  9.  Current  density  profiles  under  different  PZT  frequencies. 
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H20 

0  561 


(Rib  width  =2mm) 


Fig.  10.  Water  vapor  profiles  under  different  rib  width. 
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Fig.  11.  Oxygen  concentration  profiles  under  different  rib  width. 
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4.2.  Influence  of  rib  width  on  current  density 

Although  water  may  accumulate  near  the  rib  and  may  reduce 
internal  resistance  in  wider  ribs,  Fig.  8  shows  that  both  a  wider  rib 
and  a  thinner  rib  have  a  slight  effect  on  cell  performance.  However, 
a  current  overshoot  may  occur  at  the  initial  stage  in  the  wide  rib 
case.  The  current  will  then  decay  and  be  lower  than  in  the  thin  rib 
case. 

4.3.  Influence  of  PZT frequency  on  current  density 

As  the  PZT  vibrates  at  low  frequencies,  it  is  shown  in  Fig.  9  that 
current  generation  is  lower  compared  with  high  PZT  vibration  fre¬ 
quencies.  However,  a  current  overshoot  occurs  in  the  wide  rib  case, 
especially  at  the  high  PZT  vibration  frequency.  In  addition,  the  cur¬ 
rent  decreases  with  time  at  the  low  PZT  vibration  frequency  (/=  8, 
16  Hz)  due  to  the  lower  air-flow  rate  in  the  cathode  channel.  Fig.  10 
illustrates  that  the  current  may  decay  with  time  because  the  accu¬ 
mulated  water  vapor  cannot  be  pumped  out  immediately  in  the 
cathode  channel  with  low  PZT  vibration  frequencies.  The  simu¬ 
lations  show  that  the  water  distribution  in  the  cathode  channel 
is  only  slightly  different  under  both  rib-width  cases.  Fig.  11  also 
shows  that  the  rib  width  slightly  affects  oxygen  distribution  in  the 
cathode  channel.  Therefore,  it  may  be  the  reason  for  the  slight 
difference  in  current  density  between  the  two  different  widths 
(Fig.  9). 

The  simulation  result  of  this  study  is  a  quasi-steady  state,  and  it 
consists  of  64  cycles  in  1  s,  i.e.,/=  64  Hz.  In  order  to  reduce  compu¬ 
tational  time,  the  time  scale  is  chosen  as  1  s. 

4.4.  Influence  of  equivalence  ratio  on  current  density 

When  the  PZT  vibrates  at  high  frequencies,  more  air  can  be 
drawn  into  the  cathode  channel,  leading  to  a  reduction  in  the  equiv¬ 
alence  ratio.  Therefore,  higher  current  density  can  be  found  at 
f=  64  Hz,  where  the  equivalence  ratios  are  1.72  and  1.54  with  rib 
sizes  of  1  and  2  mm,  respectively  (Table  3).  In  contrast,  with  higher 
equivalence  ratios,  less  air  will  be  drawn  into  the  cathode  chan¬ 
nel  at  low  PZT  frequencies.  When  the  rib  width  increases  from  1  to 
2  mm,  the  difference  in  the  equivalence  ratio  is  2.76%  at  a  high  PZT 
frequency  of/=  64  Hz.  However,  the  difference  in  equivalence  ratio 
is  only  0.35%  at  a  low  PZT  frequency  of/=  8  Hz. 

The  simulations  indicate  that  the  equivalence  ratio  is  not 
affected  by  the  rib  width  (Fig.  12).  In  addition,  the  current  den¬ 
sity  is  high  at  high  frequencies.  Both  of  the  different  rib  widths 
generate  almost  the  same  current  density,  0.205  and  0.208  A  cm-2, 
using  1  and  2  mm  rib  widths,  respectively,  at  /=  64  Hz.  Therefore, 
the  effect  of  the  PZT  vibration  frequency  is  stronger  than  the  effect 
of  rib  width  in  the  air-flow  rate. 


Table  3 

Air-flow  rate  at  different  PZT  frequencies  and  equivalence  ratio  at  the  same  time 
(t=  1  s) 


Fig.  12.  Current  density  profiles  under  different  equivalent  ratio. 


5.  Conclusions 

A  novel  ribbed  PZT-PEMFCs  design  has  been  developed,  and  a 
three-dimensional,  transitional  model  has  been  successfully  built 
to  consider  its  primary  phenomena  and  fuel  cell  performance.  A 
ribbed  cathode  channel  can  reduce  the  internal  resistance  and  can 
double  the  current  density.  However,  the  ribbed  design  may  cause 
the  accumulation  of  water  vapor  near  the  rib.  At  higher  PZT  vibra¬ 
tion  frequency  (/=  64  Hz),  an  air-breathing  PZT-PEMFC  compresses 
more  oxygen  into  the  catalyst  layer  and  thus  enhances  the  electro¬ 
chemical  reaction,  resulting  in  a  higher  current  output.  At  the  same 
time,  the  accumulated  water  vapor  may  be  pumped  out  from  the 
cathode  channel  during  the  compression  process.  In  addition,  the 
effect  of  PZT  vibration  frequency  is  stronger  than  the  effect  of  the 
rib  width  on  the  air-flow  rate  and  current  density.  The  optimal  cur¬ 
rent  density  for  the  ribbed  PZT-PEMFCs  in  the  cathode  below  0.7  V 
is  0.208  A  cm-2  at/=64Hz. 
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